Structural evolution in Pt isotopes with the Interacting Boson Model Hamiltonian 
derived from the Gogny Energy Density Functional 
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Spectroscopic calculations are carried out, for the description of the shape/phase transition in 
Pt nuclei in terms of the Interacting Boson Model (IBM) Hamiltonian derived from (constrained) 
Hartree-Fock-Bogoliubov (HFB) calculations with the finite range and density dependent Gogny- 
D1S Energy Density Functional. Assuming that the many-nucleon driven dynamics of nuclear 
surface deformation can be simulated by effective bosonic degrees of freedom, the Gogny-DlS po- 
tential energy surface (PES) with quadrupole degrees of freedom is mapped onto the corresponding 
PES of the IBM. Using this mapping procedure, the parameters of the IBM Hamiltonian, relevant to 
the low-lying quadrupole collective states, are derived as functions of the number of valence nucle- 
ons. Merits of both Gogny-HFB and IBM approaches are utilized so that the spectra and the wave 
functions in the laboratory system are calculated precisely. The experimental low-lying spectra of 
both ground-state and side-band levels are well reproduced. From the systematics of the calculated 
spectra and the reduced E2 transition probabilities B(E2), the prolate-to-oblate shape/phase tran- 
sition is shown to take place quite smoothly as a function of neutron number N in the considered Pt 
isotopic chain, for which the 7-softness plays an essential role. All these spectroscopic observables 
behave consistently with the relevant PESs and the derived parameters of the IBM Hamiltonian as 
functions of N. Spectroscopic predictions are also made for those nuclei which do not have enough 
experimental E2 data. 

PACS numbers: 21.10.Re,21.60.Ev,21.60.Fw,21.60.Jz 
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I. INTRODUCTION 



The quadrupole collective motion has always attracted 
considerable attention in nuclear physics [l|-|3( . Neverthe- 
less, a fully microscopic understanding of the evolution 
of the nuclear shapes with the number of nucleons still 
remains a major challenge [3-fl3| - From the experimental 
point of view, low-lying spectroscopy is one of the most 
powerful sources of information about structural evolu- 
tion and/or shape transitions in atomic nuclei since it 
allows to establish signatures correla ting the excitation 
energies with deformation properties |1#|20| . In partic- 
ular, the complex interplay between several deformation 
degrees of freedom, taking place in different regions of the 
nuclear chart, offers the possibility of testing microscopic 
descriptions of atomic nuclei under a wide variety of con- 
ditions. In this context, mean-field approximations based 
on effective Energy Density Functionals (EDFs), which 
are a cornerstone to almost all microscopic approxima- 
tions to the nuclear many-body problem [3|, appear to 
be a first tool to rely on when looking for fingerprints of 
nuclear shape/phase transitions. 

Mean-field approximations are based on product trial 
wave functions, which are used to minimize a given EDF. 
Such products break several symmetries of the underlying 
nuclear Hamiltonian (spontaneous symmetry breaking 
mechanism) allowing the use of an enlarged Hilbert space 
within which static correlations associated with collective 



modes (e.g., quadrupole deformations) are incorporated 
at the cost of a moderate effort. Nowadays, systematic 
mean-field studies are possible because, on the one hand, 
important advances have been made in the fitting proto- 
cols providing EDFs with global predictive power all over 
the nuclear chart. Popular EDFs for calculations along 
these lines are the non-relativistic Gogny [2l|, [22[ and 
Skyrme [1, [H, [24j ones, as well as different parameteri- 
zations of the relativistic mean-field Lagrangian [f| [25[ . 
On the other hand, it has also become possible to recast 
mean-field equations in terms of efficient minimization 
procedures such as the so-called gradient method 
One of the advantages of the gradient method is the way 
it handles constraints, which is well adapted to the case 
where a large number of constraints are required (like 
the case which requires, in addition to the proton and 
neutron number constraints, constrains on both /3 and 
7 degrees of freedom characterizing the nuclear shape). 
Another advantage is its robustness in reaching a solu- 
tion, a convenient property when large scale calculations 
requiring the solution of many HFB equations are per- 
formed. 

On its own, the Interacting Boson Model (IBM) [H] 
has been quite successful in reproducing the experimen- 
tal spectra and electromagnetic transitions for low-lying 
quadrupole collective states. The virtue of the IBM is, 
with its simplicity, its robust capability of calculating 
the spectroscopic observables precisely, while the param- 
eters of the IBM Hamiltonian have been determined phe- 
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nomenologically. Therefore, the IBM itself has a certain 
microscopic foundation, where the collective J = + (S) 
and 2 + (D) pairs of valence nucleons are approximated 
by J = 0+ (s) and 2+ (d) bosons, respectively [H|. The 
proton and neutron degrees of freedom can be taken into 
account, where the so-called proton (neutron) s n and 
(.Si, and du) bosons correspond to the collective pairs 
of valence protons (neutrons) S v and D„ (S v and D v ) 
[29l [30| . This is closer to a microscopic picture compared 
to a simpler version of IBM and is known as the proton- 
neutron Interacting Boson Model (IBM-2). As the num- 
ber of valence protons (neutrons) is constant for a given 
nucleus, the number of proton (neutron) bosons, denoted 
by n-n (n v ), is set equal to half of the valence proton (neu- 
tron) numbers. The derivation of the IBM Hamiltonian 
has been studied extensively in realistic cases for nearly 
spherical or 7-unstable shapes [3lH34l using generalized 
seniority states of the shell model [29|, 30], as well as for 
deformed nuclei [35J, |36[, but still remains to be done 
for general cases in a unified manner. Therefore, it is 
timely and necessary to bridge the IBM and mean-field 
models with the help of fermion-to-boson mapping pro- 
cedures. The key question here is to investigate to which 
extent the underlying fermionic dynamics of mean-field 
models can be translated into effective bosonic degrees of 
freedom. Such an approach would enable one to take ad- 
vantage of the universality of microscopic nuclear EDFs 
d, HJ-il] and the simplicity of the IBM [Hj]. By the 
combination of both models, one would be able to access 
the spectroscopic observables which have the good quan- 
tum numbers in the laboratory system, including those 
for experimentally unexplored nuclei. 

A novel way of deriving the parameters of the IBM 
Hamiltonian has been recently proposed by two of us [13] ■ 
The IBM Hamiltonian has been constructed by mapping 
the mean- field potential energy surface (PES), obtained 
in the framework of (constrained) Skyrme Hartree-Fock 
plus BCS calculations [Hj], onto the corresponding PES 
of the IBM. The parameters of the IBM Hamiltonian, rel- 
evant to the description of the considered quadrupole col- 
lective states, have been shown to be determined uniquely 
as function of the number of valence nucleons by using 
the Wavelet analysis [39|. Calculations along these lines 
have been performed to study the shape/phase transition 
in Sm isotopes with neutron number N = 82 ~ 96, as well 
as for Ba, Xe, Ru and Pd isotopes with N — 50 ~ 82. 
Spectroscopic predictions have also been made for W and 
Os nuclei with N >126 in the lower- right quadrant of 
208 Pb [13, HI]. In addition, it has to be mentioned that 
the quantum mechanical correlation effect on the binding 
energies can be included in such calculations by diagonal- 
izing the mapped IBM Hamiltonian [39]. Note that, in 
this framework, the IBM keeps its important properties 
including the boson number counting rule as well as the 
algebraic features. 

In this paper we present, spectroscopic calculations 
for the Pt isotopic chain (i.e., for the even-even iso- 
topes 172 -20°pt) in terms of an IBM Hamiltonian de- 



termined microscopically by mapping the PES obtained 
in the framework of the (constrained) Hartree-Fock- 
Bogoliubov (HFB) approximation [t| El H3 based on 
the parametrization D1S [H of the Gogny-EDF [211221]. 
Quite recently, the structural evolution in Pt isotopes, 
including the role of triaxility (i.e., the 7 deg ree of free- 
dom), has been studied by three of us [4lj . In addi- 
tion to the (standard) Gogny-DlS EDF, the new incar- 
nations DIN and DIM gl[ of the Gogny-EDF have 
also been included in the mean- field analysis of Ref . [4l| . 
The considered range of neutron numbers included pro- 
late, triaxial, oblate and spherical shapes and served for 
a detailed comparison of the (mean-field) predictions of 
the new parameter sets DIN and DIM against the stan- 
dard parametrization D1S. It has been shown that, re- 
gardless of the particular version of the Gogny-EDF em- 
ployed, the prolate-to-oblate shape/phase transition oc- 
curs quite smoothly with the 7-softness playing an im- 
portant role. It is therefore very interesting to study how 
the systematics of the HFB PESs discussed in Ref. [4l[ 
is reflected in the isotopic evolution of the corresponding 
low-lying quadrupole collective states and how accurately 
such states can be reproduced by a mapped IBM Hamil- 
tonian [13, Hi) ■ Let us stress that our main goal in the 
present work is to study the performance of a fermion-to- 
boson mapping procedure [37], Hi] based on the Gogny- 
EDF. For this reason, as a first step, we will restrict 
ourselves to a mapping in terms of the parametrization 
Gogny-DlS already considered as global and able to de- 
scribe reasonably well low-energy experimental data all 
over the nuclear chart (see, for example, Refs. [13, |H| 
and references therein). 

From the theoretical perspective, the Pt and neighbor- 
ing isotopic chains have been extensively studied in terms 
of both IBM and mean-field-based approaches. There is 
much experimental evidence [45l |46| revealing existences 
of 7-unstable 0(6) nuclei in Pt isotopes. The IBM-2 
has been used in a phenomenological way for the spec- 
troscopy of Pt, Os and W isotopes [13, S3- The prolate- 
to-oblate transition in Pt as well as in Os and W nuclei, 
has been observed in the recent experiment [49| , where a 
relatively moderate oblate-to-prolate shape/phase transi- 
tion occurs in Pt as compared to Os and W nuclei. Spec- 
troscopic calculations have been carried out for Pt iso- 
topes in the framework of the five-dimensional collective 
Hamiltonian, derived from the pairing-plus-quadrupole 
model 50]. Evidence for 7 vibrations and shape evo- 
lution in 184 - 190 Hg has been considered in Ref. [5lj ]. 
where a five-dimensional collective Hamiltonian was built 
with the help of constrained Gogny-DlS HFB calcula- 
tions. On the other hand, systematic mean-field studies 
of the evolution of the ground state shapes in Pt and the 
neighboring Yb, Hf, W and Os nuclei have been carried 
out with non-relativistic Skyrme [1] and Gogny [23, El| 
EDFs, as well as within the framework of the relativis- 
tic mean- field (RMF) approximation [52|. One should 
also keep in mind, that Pt, Pb and Hg nuclei belong to 
a region of the nuclear chart, around the proton shell 
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closure Z — 82, characterized by a pronounced compe- 
tition between low-lying configurations corresponding to 
different intrinsic deformations [H3| and therefore, a de- 
tailed description of the very rich structural evolution in 
these nuclei requires the inclusion of correlations beyond 
the static mean-field picture [H3 - l56l | accounting for both 
symmetry restoration and configuration mixing. The role 
of configuration mixing in this region has also been con- 
sidered in phenomenological IBM studies [13, HI] . 

The paper is organized as follows. In Sec. HH we will 
briefly describe the theoretical tools used in the present 
study. Illustrative examples of IBM PESs, obtained by 
mapping the corresponding Gogny-HFB PESs, are pre- 
sented in Sec. IIIII The isotopic evolution of the IBM 
parameters derived for the nuclei 172 - 200 pt is discussed 
in Sec. IIV1 Spectroscopic calculations, including the sys- 
tematics of excitation spectra and reduced E2 transition 
probabilities B(E2) along the Pt isotopic chain, will be 
discussed in Sec. |V] There, we will also show detailed 
comparisons between the predicted level schemes and the 
available data for some Pt isotopes selected as a repre- 
sentative sample. Finally, Sec.|VT]is devoted to the con- 
cluding remarks and work perspectives. 



II. THEORETICAL PROCEDURE 

In this section, we briefly describe the theoretical 
frameworks used in the present study, i.e., the con- 
strained HFB approximation, as well as the procedure 
followed to construct the corresponding mapped IBM 
Hamiltonian. For more details the reader is referred to 
Refs. [Him and [H. 

In order to compute the Gogny-HFB PESs, which are 
our starting point, we have used the (constrained) HFB 
method together with the parametrization D1S of the 
Gogny-EDF. The solution of the HFB equations, leading 
to the set of vacua |$hfb(/3, 7)), is based on the equiv- 
alence of the HFB with a minimizati on p roblem that is 
solved using the gradient method [26|, [27|. In agreement 
with the fitting protocol of the force, the kinetic energy of 
the center of mass motion has been subtracted from the 
Routhian to be minimized in order to ensure that the 
center of mass is kept at rest. The exchange Coulomb 
energy is considered in the Slater approximation and we 
neglect the contribution of the Coulomb interaction to 
the pairing field. The HFB quasiparticle operators are 
expanded in a Harmonic Oscillator (HO) basis contain- 
ing enough number of shells (i.e., N s h e ii = 13 major 
shells) to grant convergence for all values of the mass 
quadrupole operators and for all the nuclei studied. We 
constrain the average values of the mass quadrupole op- 
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to the desired deformation values Q20 ancl Q22 defined as 



Q20 = (3>hfb|<32o|3>hfb) 
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Q22 = ($hfb|Q22| < I , hfb)- (2) 
In Ref . [4l[ , the Q — 7 energy contour plots with 
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have been used to study the (mean-field) evolution of 
the ground state shapes in Pt nuclei. Alternatively, one 
could also consider the (3 — 7 representation in which the 
quadrupole deformation parameter /3 is written [27| in 
terms of Q [Eq.©] as 



P = 



Q 



A(r 2 ) 



(5) 



where (r 2 ) represents the mean squared radius evaluated 
with the corresponding HFB state |$hfb)- 

The set of constrained HFB calculations described 
above, provides the Gogny-DlS /3 — 7 PES (i.e., the to- 
tal HFB energies £hfb(/#, 7) 0) required for the subse- 
quent mapping procedure, for which the following IBM-2 
Hamiltonian Hibm is employed 



where 



and 



ndp 



= <#-d p , (p = TT,v) 



Qp = [std p + d\s p ] {2) + Xp^pdpi 



1(2) 



(6) 



(7) 



(8) 



stand for the <i-boson number operator and the 
quadrupole operator, respectively. The competition be- 
tween the coupling constants e and n determines the de- 
gree of nuclear deformation. 

The bosonic PES is represented by the expectation 
value of Hibm, computed in terms of the so-called bo- 
son coherent state [59M6lj | 



n [* 



p— 7T,I/ 



p=0,±2 



|o> 



(9) 



and 



where |0) stands for the boson vacuum (i.e., inert core) 
and the coefficients a's are expressed as a p o = /3 p cos7 P , 
a p ±i — and a p ±2 — -^/? p sin7 p . Within this context, 
the intrinsic shape of the nucleus is described in terms of 



4 



the (axially symmetric) deformation j3 p and the (triaxial) 
deformation j p . In the present study, as well as in our 
previous works [37], l39j |. we assume for simplicity that 
= f3 u =J B and 7^ = j v = 7 B - The IBM PES is then 
given by HEH 



^ibm(/?b,7b) — 



1 + 



(1 + /3I) 2 



4 - 2\ -(xtt + X»)Pb cos 37b + -=X-kXv0b 



(10) 



Here we assume the proportionality /?b = C/3/3, with 
(7g being a numerical coefficient [37]. If one further as- 
sumes the separability of the mapping along the /3 and 7 
directions [13, HH, one then has 7b = 7- Thus, (/3b, 7b) 
represent the boson images of the (fermion) deformation 
parameters ((3, 7) given by Eqs. (@| and ([5]). We then 
map a point on the HFB PES, (,9,7), within an energy 
range relevant for the considered low-lying quadrupole 
collective states, onto the corresponding point on the 
IBM PES, (/9b ,7b)- This process is exactly the mapping 
of the fermionic PES onto the bosonic one. In practice, 
one determines the 6, ft, Xtt,^ 

and Cp values for each in- 
dividual nucleus by drawing the IBM PES so that the 
topology of the corresponding HFB PES is reproduced. 
This is done unambiguously by means of the recently de- 
veloped procedure [39|, which makes use of the powerful 
method of the Wavelet transform [62| . 

Here we would like to make the following remarks: The 
topology of the HFB PES reflects essential (fermionic) 
features of many-nucleon systems, such as the Pauli prin- 
ciple and the underlying nuclear interactions. Such ef- 
fects are supposed to be incorporated into the boson sys- 
tem by the mapping procedure [37J. On the other hand, 
the solution of the five-dimensional (5D) collective Bohr 
Hamiltonian with parameters obtained from EDFs cal- 
culations (see, for example, Q [H, [641 ]) is a popular 
alternative to obtain the low-lying collective spectra in 
even-even nuclei. In this kind of calculations the pure 
mean-field PES is replaced by another quantity that in- 
corporates in addition to the HFB energy, the zero point 
rotational and vibrational corrections. These corrections 
to the energy are intimately related to the use of a gener- 
alized kinetic energy term for the collective motion that 
includes not only moments of inertia but also collective 
vibrational masses. To what extent the present mapping 
procedure plus the solution of the IBM Hamiltonian is 
able to mimic the solution of the 5D Bohr Hamiltonian 
is still an open question, that can be partially answered 
by looking at the reasonable results obtained with our 
method and that compare qualitatively well with the ones 
of the 5D Hamiltonian. A possible way to incorporate the 
effect of the collective masses into the mapping would be 
to make a change of variables (analogous to the one in- 
voked in the derivation of the GOA Q) as to render the 
collective masses constant all over the range of allowed 
values of the j3 and 7 deformation parameters. Perhaps, 
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this is the missing element that could correct some ob- 
served [37|, l39| systematic deviations of the IBM rota- 
tional spectra with respect to the experimental ones for 
well deformed systems, and requires the introduction of 
an additional mass term in the IBM Hamiltonian known 
as the L ■ L term (65|. However, this problem does not 
show up for moderately deformed cases like the ones stud- 
ied in the present work and therefore it is not considered 
here. 



III. MAPPED POTENTIAL ENERGY 
SURFACES 

The IBM PESs obtained for the nuclei 1 80-i98p t are 
shown in Fig. Q] as a representative sample. The IBM 
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parameters e, k, Xir,v and Cp, to be discussed later on in 
Sec. IIV[ have been obtained by mapping the correspond- 
ing Gogny-DlS PESs presented in Fig. 2 of Ref. [Zl| along 
the lines previously described in Sec. [TTJ 

The IBM PESs from 180 Pt to 186 Pt display a prolate 
deformed minimum and an oblate deformed saddle point. 
The prolate minimum becomes softer in 7 but steeper in 
/3 direction as the number of neutrons increases. This is, 
roughly speaking, consistent with the topologies of the 
HFB PESs of Ref. [4l| , where the minima are located a 
bit off but quite nearby the line 7 = 0°. 

The IBM PESs for both 1 88,i90p t are r) soft; having the 
minimum on the oblate side. These nuclei are supposed 
to be close to the critical point of the prolate-to-oblate 
shape transition. The corresponding HFB PESs display 
shallow triaxial minima with 7 ~ 30° and are also soft 
along the 7 direction [4l|. The IBM Hamiltonian con- 
sidered in the present study does not provide a triaxial 
minimum, but either prolate or oblate minimum, as can 
be seen from Eq. (TTU)) . The 7-softness can be simulated 
by choosing the parameters Xir and Xv so that their sum 
becomes nearly equal to zero. This is reasonable when 
a triaxial minimum is not deep enough like the present 
case, where the triaxial minimum point in the HFB PES 
differs by at most several hundred keV in energy from 
either prolate or oblate saddle point. However, the topol- 
ogy of the mapped IBM PES is then somewhat sensitive 
to the values of the parameters Xn and Xvi which occa- 
sionally results in a quantitative difference in the loca- 
tion of the minimum in the IBM PES from that of the 
HFB PES. In fact, and contrary to what happens with 
the HFB PESs plj], the IBM PES of 190 Pt is softer in 
7 than that of 188 Pt. One should then expect a certain 
deviation of the resultant IBM spectra from the experi- 
mental ones, which can be partly attributed to the small 
difference already mentioned. 

In Fig. [TJ isotopes from 192 Pt to 198 Pt exhibit oblate 
deformation. The locations of their energy minima and 
their curvatures in both j3 and 7 directions agree well 
with the ones of the Gogny-DlS PESs in Ref. HlJ. These 
isotopes become steeper in the 7 direction and shallower 
in the /3 direction as the number of neutrons increases. 
Their energy minima approach the origin more rapidly 
than the lighter Pt nuclei shown Fig. [TJ This evolution 
reflects the transition from oblate deformed ground states 
to a spherical vibrator as one approaches the neutron 
shell closure N = 126. 



IV. DERIVED IBM PARAMETERS 



The IBM parameters e, n, Xn,v and Cp derived for the 



nuclei 



Pt from the mapping procedure described 



in Sec. [H] are depicted in Figs. [Ua)-[21]d) as functions of 
the mass number A. 

Figure [Ha) shows the parameter e gradually decreases 
toward mid shell in accordance with the growth of the 
deformation. This trend reflects the structural evolution 
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FIG. 2: (Color online) IBM parameters e, k, X^<" an d C73 as 
functions of the mass number A. For the wavelet analysis the 
Morlet function is used [391 ■ 



from nearly spherical to more deformed shapes and is 
consistent with previous results for other isotopic chains 
[3!^ ]. In Fig. djb), the derived k parameter is almost 
constant and somewhat larger in comparison with the 
phenomenological value [48| , which is the consequence of 
the sharp potential valleys observed in the Gogny-DlS 

PESs jiy. 

On the other hand, in Fig. [21c) the proton parame- 
ter Xir is almost constant while the neutron parameter 
Xv changes significantly. The systematic behavior of the 
present Xv value is consistent with the phenomenological 
one [48j], while there is quantitative difference between 
the former and the latter. The magnitude of the sum 
Xn + Xv as well as its sign depend on how sharp the HFB 
PES is in the 7 direction and on whether the nucleus is 
prolate (negative sum) or oblate (positive sum) deformed, 
respectively. Therefore, as Xir does not change much, the 
role of 7 instability can be seen clearly from the system- 
atics of Xv For the isotopes 172 - 180 Pt the PES exhibits 
prolate deformation and the sum x-n + Xv has negative 
sign. The average of the derived x-rr and Xv values is 
nearly equal to zero for the nuclei 182 - 194 pt. This is a 
consequence of the 7 softness in the corresponding HFB 
PESs. On the other hand, the sum Xtt + Xv becomes 
larger with positive sign as we approach the neutron shell 
closure iV=126 reflecting the appearance of weakly de- 
formed oblate structures in the corresponding PESs. 

Figure [^d) shows that Cp decreases gradually toward 
the middle of the major shell. Cp can be interpreted as 
the "bridge" between the geometrical deformation /3 [l[ 
and the IBM deformation /3b and is thus proportional to 
the ratio between the total and valence nucleon numbers, 
in a good approximation [60] . This is probably the reason 
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for the decreasing trend observed in Fig.[2fd), as well as 
in earlier studies for other isotopic chains[37], |39| . 



V. SPECTROSCOPIC CALCULATIONS 

With all the parameters e, k, x-k,v and Cp required 
by the IBM Hamiltonian at hand, we are now able to 
test the spectroscopic quality of our mapping procedure, 
based on the Gogny-DlS EDF, for the nuclei i72-200p t _ 
Therefore, in the following we will discuss our predic- 
tions concerning the properties of the low-lying spectra 
as well as the reduced transition probabilities B(E2). We 
will also consider their correspondence with the mapped 
PESs and the derived IBM parameters. We will compare 
our theoretical predictions with the available experimen- 
tal data taken from Brookhaven National Nuclear Data 
Center (NNDC) [66j and from the latest Nuclear Data 
Sheets [67(. The diagonalization of the IBM Hamilto- 
nian is performed numerically for each nucleus using the 
code NPBOS [H. 

Here we have to note that the experimental 2g~ and 
4^" levels for mass numbers 192 ^ A ^ 200 belong to 
bands different from that of the 0^" level, while they are 
assigned to the quasi-/3-band levels in Fig. IS^b), as well 
as in Fig. [SJ for convenience sake. Similarly, as one will 
see in Fig. [3tc) the experimental data for the 3^ and the 
Pt are assigned to the quasi-7-band 



4+ levels in 198 - 200 



levels lying on top of the 2 J energy. 



A. Evolution of low-lying spectra 

Figure |3] displays the calculated spectra for (a) ground- 
state, (b) quasi-/? and (c) quasi-7 bands. What is striking 
is the good agreement between the present calculations 
and the experimental data not only for ground-state but 
also for quasi-/3 and quasi-7 band energies, where overall 
experimental trends are reproduced fairly well in partic- 
ular for the open-shell nuclei 18 °- 192 pt. 

We show in Fig. [3Ja) the evolution of the 2f, if, 6+ 
and 8^ levels in the considered Pt nuclei as functions of 
the mass number A. The calculated energies decrease 
toward the middle of the major shell with the number 
of the valence neutrons and remain almost constant for 
176 < A < 186 nuclei. Although these tendencies are 
well reproduced, the rotational features are somewhat 
enhanced in the calculated levels for 180 > 182 > 184 pt which 
are slightly lower in energy than the experimental ones. 
From both the theoretical results and the experimental 
data, one can also observe clear fingerprints for structural 
evolution with a jump between 186 Pt and 188 Pt, which 
can be correlated with the change of the mapped PESs 
from prolate to oblate deformations. For A ^ 188 the 
yrast levels gradually go up as the neutron shell closure 
iV=126 is approached. 

One can also find signatures for a shape/phase transi- 
tion in the systematics of the quasi-/? band levels shown 



in Fig.EJb). From A =180 to 186, the 0^ band head and 
the 2^" level look either constant or nearly constant in 
both theory and experiment. The two levels are pushed 
up rather significantly from A=186 to 188 consistently 
with the systematics in the ground-state band and with 
the change of the mapped PESs as functions of the neu- 
tron number N. The calculated 0^~ and 2^" levels are 
higher than but still follow the experimental trends. 

Coming now to the quasi-7 band levels shown in 
Fig. Etc) , one can observe the remarkable agreement be- 
tween theoretical and experimental spectra for 180 ^ 
A ^ 186, where the 3+ level lies close to the 4^ level. 
However, the present calculation suggests this trend per- 
sists even for 188 ^ A ^ 196, whereas the relative spacing 
between the experimental and 4^" levels for these nu- 
clei is larger. Similar deviation occurs for 5^ and 62" lev- 
els, although the latter is not exhibited in Fig.^b). This 
means that our calculations suggest the feature charac- 
teristic of the 0(6) symmetry, where the staggering oc- 
curs as 2+, (3+ 4+), (5+ 6+), .... etc. However, the 
experimental levels are lying more regularly particularly 
for 188 ^ A ^ 196, and thus appear to be in between 
the 0(6) limit and a rigid triaxial rotor where the stag- 
gering shows up as (2+ 3+), (4+ 5+), ... etc [69]. Such 
a deviation of the 7-band structure seems to be nothing 
but a consequence of an algebraic nature of the IBM, and 
indeed has also been found in existing phenomenological 
IBM calculations [48[ . From a phenomenological point of 
view, the so-called cubic (or the three-boson) interaction 
[7Cl l71j has been useful for reproducing the experimental 
7-band structure. The cubic term produces a shallow tri- 
axial minimum that is seen in the Gogny-HFB PES, and 
may be introduced also in the Hamiltonian in Eq. ([5]). 
This is, however, out of focus in the current theoretical 
framework, because the cubic term represents an effective 
force whose origin remains to be investigated further. 

Here, the deviations observed in the side-band lev- 
els (even in some of the ground-state band levels) for 
A 196 are probably related to the larger magnitude of 
the parameter k as compared with its phenomenological 
value [481 ] . Roughly speaking, when the magnitude of k 
becomes larger, the moment of inertia decreases, result- 
ing in the deviation of not only ground-state-band but 
also the side-band energies. The problem arises in the 
present case partly because, in the vicinity of the shell 
closure N = 126, the HFB PESs exhibit weak oblate de- 
formations close to the origin j3 = 0, as we showed in 
Fig. [1] In addition, the curvatures along the j3 direction 
around the minima are somewhat larger. These peculiar 
topologies of the Gogny-DlS PESs make it rather dif- 
ficult to determine a value of k which gives reasonable 
agreement of side-band energies with the experimental 
ones. In this case one may interpret that the deviation 
is mainly due to the properties of the particular version 
of the Gogny-EDF considered in the present study. An- 
other possibility is that the boson Hamiltonian used may 
be still simple, requiring the introduction of additional in- 
teraction terms in the boson system. Investigation along 
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FIG. 3: (Color online) Evolution of calculated (curves) and experimental (symbols) low-lying spectra of 172 200 Pt nuclei for 
(a) ground-state, (b) quasi-/? and (c) quasi-7 bands as functions of the mass number A. Experimental data are taken from 
Ref. [Hi] 



these lines is in progress and will be reported elsewhere. 



B. Systematics of J3(E2) ratios 

Once the boson wave functions corresponding to the 
excited states of a given nucleus are obtained, we are able 
to compute electromagnetic transitions, among which the 
reduced E2 transition probabilities B(E2) are of partic- 
ular importance. The B(E2) transition probabilities are 
given by pj 



B(E2; J -> J') = — -\{J'\\f < E2 )||J)| 2 , (11) 



where J and J' are the angular momenta for the initial 
and final states, respectively. The E2 transition opera- 
tor f ( E2 > is given by f ( E2 ) = e^Q* + e v Q v , with and 
e v being the boson effective charges. In principle, the 
effective charges should be determined independently of 
the underlying mean-field calculation. In the present pa- 
per, we assume — e„, for simplicity, and focus our 
discussion on the B(E2) ratios defined as 



Ri 

R2 
R 3 



B(E2;4{ 
B(E2;2+ 
S(E2;0+ 
B(E2; 2+ 



2+)/B(E2;2+ 
2+)/B(E2;2+ 
2+)/B(E2;2+ 
0+)/B(E2;2+ 



2t 



(12) 



which arc shown in Fig. [4] as functions of the mass number 
A. 

The ratio R\ is nearly constant all the way, being much 
below the U(5) limit of IBM (Ri = 2), and is rather close 
to i?i = 10/7, which is the 0(6) and SU(3) limit of IBM. 
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FIG. 4: (Color online) Reduced E2 transition probabilities 
B(E2) of 172 - 200 pt nuclei as functions of the mass num- 
ber A. (a) _B(E2) values for low-lying states normalized by 
B(E2;2+ -> 0+) value, R u R 2 and R 3 . (b) B(E2) branch- 
ing ratio, R 4 = B(E2;2+ ->■ 0+)/B(E2;2+ -> 2+). Ratios 
Ri — R4 are defined in Eq. (|12[) . Experimental data are taken 
from Ref. m7|]. For more details, see the main text. 



Thus, Ri is not a sensitive observable to distinguish be- 
tween axially symmetric and 7-soft nuclei. This is rea- 
sonable because the structural evolution between axially 
symmetric deformed and the 7-unstable shapes is shown 
to take place quite smoothly from the systematics of the 
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mapped PESs (in Sec. and the derived IBM param- 
eters (in Sec. IIV[) . The flat behavior of Ri value for Pt 
isotopes in Fig. |4] differs from the one found e.g., in Sm 
isotopes There, a sharp decrease of R\ value can be 
seen in the line of U(5)-SU(3) shape/phase transition. 

One can see that, in contrast to the flat systematics 
of the R\ value with respect to the mass number A, the 
ratio i?2 changes significantly and is relatively large for 
i86-i96p t nuc i eii being close to ^ (0(6) limit). This is 
consistent with the softness of the PESs for these nu- 
clei. Therefore, the quantity R2 is quite sensitive to 
the shape evolution encountered in the PESs and can 
be thus considered as the best signature for 7-softness 
among R1-R4. There are not much available data over- 
all, but the experimental Ri value is also relatively large 
around 192 Pt. For the nuclei 176 - 184 pt ; the theoretical 
i?2 value is close to zero (the SU(3) limit) and slightly 
goes up from A=174 to 172, probably approaching the 
U(5) vibrational limit (i?2=2) in the vicinity of the neu- 
tron shell closure N — 82. 

Unlike the R2 case, the calculated ratio R3 does not 
change much with mass number A and is close to zero 
(0(6) and SU(3) limits of R 3 ) for 188-I96p t From 
A=180, the i?3 turns to increase as we move towards the 
neutron shell closure N = 82 and is expected to approach 
the U(5) limit (R3 = 2). The calculated R3 value is, how- 
ever, still much smaller than the experimental value at 
A = 198. In fact, both the HFB and the mapped PESs 
for the nucleus 198 Pt display a weakly deformed shape, 
which somewhat differs from the vibrational feature ex- 
pected from the corresponding experimental levels. The 
present R3 value does not exhibit a drastic change ob- 
served in shape transitions in A ^130 Ba-Xe and A ^100 
Ru-Pd isotopes, where the E2 transition from the Of 
state to the 2f is much enhanced [39|. 

Finally, the branching ratio R4 also corresponds to a 
gradual shape transition. The present calculations sug- 
gest that the R4 value is nearly zero (0(6) limit) in the 
region where the nuclei are soft and where the R2 ratio 
takes large values. The calculated R4 ratio becomes rela- 
tively larger for A ^ 184, where the PESs show stronger 
prolate deformation. Consistently with the evolution of 
the IBM PESs, the calculated R4 values turn to approach 
the U(5) limit, which is also zero, for A ^ 178. Similarly 
to the R3 deviation from the vibrational character 

of the experimental data is found at A = 198. 

It should be emphasized that all the results for _B(E2) 
values shown so far are quite consistent with the topolo- 
gies of the PESs and with the derived IBM parameter 
values. 



C. Level schemes of selected nuclei 

As already mentioned in Sec. U one of the main goals 
of the present study is to test the spectroscopic qual- 
ity of the mapping procedure and the underlying (uni- 
versal) Gogny-DlS EDF [13, EJ], which have been de- 



scribed in Sec. ITT] Keeping this in mind, we will now 
turn our attention to a more detailed comparison be- 
tween our results and the available experimental data 
for excitation spectra and B(E2) values. To this end, 
we select the nuclei i84-i94 pt 

as a representative sample 
corresponding to the mapped PESs shown in Fig.[T] The 
level schemes obtained for the nuclei 184 -!94pt are com- 
pared in Fig. [5] with available experimental data. The 
theoretical B(E2) values are shown also in Fig. [5] Note 
that B(E2;2f — > Of) value is normalized to the exper- 
imental one. The virtue of the present calculation is to 
give predicted B(E2) values for those nuclei which have 
no enough E2 information available. This is particularly 
useful in the cases of 184 Pt, 186 Pt and 188 Pt where the 
calculated spectra agree well with the experiment. 

For clarity, we divide the explanation of the results 
shown in Fig. [5] into the following three categories, ac- 
cording to the tendencies of what we found in the PESs 
and the IBM parameters and the experimental data. 
The first is the prolate deformed regime represented by 
the nuclei 184 > 186 Pt, which exhibit a rotational charac- 
ter. Next, we will consider the isotopes 188 > 190 Pt which 
are apparently close to the critical point of the prolatc- 
to-oblate transition observed in the mapped PESs (see 
Fig. [I} . Lastly, calculated and experimental results are 
compared for the nuclei 192 ^ 194 pt which belong to the 
weakly oblate deformed regime. Note that, in Fig.[5j the 
energy scale is not common for all nuclei. 

For 184 ,i 8 6p^ t ne present calculation reproduces over- 
all pattern of the experimental spectra in all of the 
ground-state, quasi-/? and quasi-7 bands fairly well. In- 
teresting enough, the bandhead energies, particularly 
the quasi- j3 bandhead Of , are much higher than the Af 
level, compared to the experimental data. This indicates 
that, reflecting the topologies of the PESs in Fig. [TJ the 
18 ,186 Pt nuclei deviate from the 7-soft 0(6) character 
and exhibit rather rotational features. The Of energies 
are predicted to be above 4^ ones in both nuclei. 

On the other hand, both 188 Pt and 190 Pt, whose PESs 
are quite flat along the 7 direction in Fig. [TJ appear to 
be closer to the 7-unstable 0(6) limit of the IBM than 
the two nuclei already mentioned above. The 2f and 4f 
levels lie close to each other in the present study and, in 
the spirit of group theory, are supposed to have the same 
t — 2 quantum number of the 0(6) dynamical symmetry 
[28| . Similarly, in our calculations the 6^, 4f , 3^ and 
Of levels are almost degenerate and can be then grouped 
into the r = 3 multiple. Along these lines, we can observe 
characteristic E2 decay patterns that are quite consistent 
with the At = ±1 selection rule of the 0(6) limit }28j . 
For instance, the transition from the Of level (supposed 
to have r = 3) to 2f level (supposed to have r = 2) is 
dominant over the one to 2f (supposed to have r = 1) in 
both 188 Pt and 190 Pt. The trend characteristic of 0(6) 
symmetry is clearly seen particularly in 190 Pt, where the 
sum of the parameters x-k and \ v almost vanishes as seen 
from Fig.[5Jc). This means that the nucleus is close to the 
pure 0(6) limit, and is consistent with the mapped PES 
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FIG. 5: (Color online) Level schemes for 1S4 ~ 194 pt nuclei. The Gogny D1S interaction is used. The theoretical B(E2) values 
(in Weisskopf units) for each nucleus are normalized to the experimental B(E2;2f — > 0±) value. Note that, in 190 Pt, the 
theoretical 3^ level is almost identical to but is lower by only 10 keV in energy than the 4^" level, where the E2 transitions 
shown are those from the 3^ state, not from the 4^ state. For details, see the main text. 



in Fig. [T] that is nearly flat along the 7 direction. Never- 
theless, the structure of the corresponding experimental 
7 band appears to have a more triaxial nature, where the 
3^ and the 4^~ levels are apart from each other. As we 
have anticipated in Sec. IIII1 this deviation arises partly 
due to the difference of the position of energy minimum 
between Gogny-DlS PESs of [4l| and the corresponding 
IBM PESs of Fig. CD 

For 192 ' 194 pt nuclei, the theoretical 7-band structure 
still looks like that of 0(6) symmetry What is of particu- 
lar interest here is that, for both 192,194 Pt, the relative lo- 
cation of the quasi-/3-band head 0^" energy is reproduced 
fairly well lying close to the 4^~ level. In addition, for 
194 Pt, the present calculation suggests that the 0^ — > 2\ 
E2 transition is dominant over the 0^ ->■ 2+ E2 tr an- 
sition, which, although there is quantitative deviation, 
agrees with the experimental trend. The reason why such 
a quantitative difference occurs may be discussed in the 
future. Compared to the experimental data, the theoret- 
ical quasi-7 band is rather stretched and the band head 
2j energy is somewhat large. The calculated 0^" energy 
is also higher than the experimental one in particular 
for 192 Pt. Accordingly, the theoretical B(E2;2^ -)■ 2\) 
value is much smaller than experimental value with re- 
spect to the -B(E2; 2^ — > Of) value. The deviations occur 
due to the derived n value, which is somewhat larger than 
the phenomenological one [48|. For relatively high- lying 



side-band 2^~ and 4^" energies, the calculated results may 
not seem to be much reliable, because even the ordering 
of these levels are not reproduced for 192 Pt. 

VI. SUMMARY 

To summarize, spectroscopic calculations have been 
carried out, for the Pt isotopic chain in terms of the Inter- 
acting Boson Model Hamiltonian derived microscopically 
based on the (constrained) Hartree-Fock-Bogoliubov ap- 
proach with the Gogny-DlS Energy Density Functional. 

The Gogny-HFB calculations provide the potential 
energy surface (PES), which reflects, to a good ex- 
tent, many-nucleon dynamics of surface deformation with 
quadrupole degrees of freedom and structural evolution 
in a given isotopic chain. By following the procedure 
proposed in Ref. |3, the PES of the Gogny-DlS EDF is 
mapped onto the corresponding bosonic PES, and can be 
then utilized as a guideline for determining the param- 
eters of the IBM Hamiltonian. This enables one to cal- 
culate the spectroscopic observables with good quantum 
numbers (i.e., the angular momentum and the particle 
number) in the laboratory system without adjustment of 
levels. 

By this approach, global tendencies of the experimen- 
tal low-lying spectra of 172 - 200 Pt nuclei are reproduced 
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quite well not only for ground-state but also for side 
bands of mainly open-shell nuclei. It has been shown 
that shape/phase transition occurs quite smoothly from 
prolate to oblate deformations as a function of N in the 
considered nuclei 186 - 192 pt ) where the 7 instability plays 
an essential role. From the analysis in Fig. [TJ the change 
of the mapped IBM PESs in 7 direction has been more 
vividly seen than in j3 direction, similarly to the corre- 
sponding Gogny-HFB PESs of Ref. [4l[ . This is consis- 
tent with the conclusions in our earlier work [4l| and also 
with many others along the same line. We have shown 
that the calculated spectra and the B(E2) ratios behave 
consistently with the evolution of the topologies of the 
mapped PES's and with the systematics of the derived 
IBM parameters as functions of the neutron number N. 
These derived parameters are qualitatively quite similar 
to the existing phenomenological IBM studies 0, l48j . 
By studying the level schemes in detail in comparison 
with the available experimental data, the present calcu- 
lation agrees with the data fairly nicely and reflects the 
algebraic aspects of the IBM, e.g., the At = ±1 selection 
rule of the E2 decay patterns. We have also made pre- 
dictions on some E2 transition patterns. These behaviors 
of the £?(E2) may need to be examined experimentally 
particularly for lighter, A < 190 nuclei with which there 
is currently few available data. 

The evolution of ground-state shape as a function of 
both N and Z has been studied within neig hboring iso- 



topic chains such as Os, W, Hf and Yb [27[. More sys- 
tematic analysis is in order for these nuclei, by more ex- 
tensive application of the present approach. It should be 
then of interest to study how the corresponding spectra 
and transition probabilities behave. 



On the other hand, the IBM Hamiltonian of Eq. ([6]) 
has rather simple form consisting of single-d-boson oper- 
ator and the quadrupole-quadrupole interaction between 
proton and neutron bosons. The results provided by 
the present Hamiltonian were shown to be already quite 
promising. However, more studies may be necessary in 
the future for further refinement, e.g., in describing de- 
tailed structure of the quasi-7 band. Work along this line 
is in progress. 
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